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A density functional theory of ferromagnetism in heterostructures of compound semiconductors 
doped with magnetic impurities is presented. The variable functions in the density functional 
theory are the charge and spin densities of the itinerant carriers and the charge and localized spins 
of the impurities. The theory is applied to study the Curie temperature of planar heterostructures of 
III-V semiconductors doped with manganese atoms. The mean-field, virtual-crystal and effective- 
mass approximations are adopted to calculate the electronic structure, including the spin-orbit 
interaction, and the magnetic susceptibilities, leading to the Curie temperature. By means of these 
results, we attempt to understand the observed dependence of the Curie temperature of planar 5- 
doped ferromagnetic structures on variation of their properties. We predict a large increase of the 
Curie Temperature by additional confinement of the holes in a 5-doped layer of Mn by a quantum 
well. 



PACS numbers: 75.70.Cn, 75.50.Pp, 75.10.-b 
I. INTRODUCTION 

Interest in ferromagnetic III-V semiconductors lies 
both in fundamental physics and potentially useful tech- 
nological applications utilizing spins.u The growth. 
Gai-^Mn^As, a ferromagnetic III-V semiconductor,! 
has raised the basic problems of the origin of the ferro- 
magnetism and of the spin transport properties. In coru 
trast to the much studied Mn-doped II- VI materialsjD 
Mn acts as an acceptor in GaAs so that Gai-^Mn^As 
has free holes which are thought to be responsible for 
the high Curie temperaturea (Tc) of 110 K for x = 
0.054. Progress in "spintronics" is made by the recent 
demonstrations of the injection of a spin-polarized cur- 
rent from both, ferr omagnetic metalsBllj and magnetic 
semiconductorst3li3 into a semiconductor. 

The Zener model of fer roma gnetism for a bulk alloy 
of (III-V) semiconductors Ji3li3 linking the Curie Tem- 
perature with the spin susceptibility of the mobile holes, 
forms a framework for understanding the dependence of 
the ferromagnetism on various properties of the system, 
including the hole density and the Mn concentration. 
The semiconductors afford an opportunity to change the 
properties to affect the ferromagnetism, for example, 
by optical excitation of the carrierst^l and by field ef- 
fect control. O Engineering the band gap profiles of pla- 
nar heterostructures of the III-V semiconductors and the 
doping can also vary the factors influencing the ferromag- 
netic order. In experiments^ in a double quantum well 
of GaAsMn/GaAsAl, the coupling between the two mag- 
netic layers is observed to depend on both the thickness 
and the composition of the non-magnetic barrier. 

A different kind of planar system, the so-called digital 
ferromagnetic heterxistructure, has been introduced by 
R. Kawakami et alx3 This system consists of a sequence 
of atomic monolayers of Gai-^Mn^As (with x=0.25 and. 
0.5) separated by several layers of GaAs. This (5-dopingE2l 
structure displays a transition temperature which is a de- 



creasing function of the distance between the magnetic 
layers— At the shortest distance reported, Tc is around 
50 K.I13 At large interlayer distance, the layers are decou- 
pled and Tc reaches a value of arouad 35 K. Single layer 
samples also show ferromagnetism.E2l 

In this paper we provide a theoretical framework to 
study ferromagnetism in planar heterostructures of ferro- 
magnetic III-V semiconductors. The model combines the 
standard procedure to calculate the electronic structure 
of planar semiconductor heterostructures in the envelope 
function formalism with a mean field theory for the ferro- 
magnetic state. Our goal is to understand the interplay of 
confinement, spin-orbit interaction, and close packing of 
Mn atoms in the ferromagnetic digital heterostructures. 
Keeping the theories of the electronic structure and of 
ferromagnetism simple enable us to study the effects of 
varying the system configurations, degree of interdiffu- 
sion, and carrier compensation. We apply the theory to 
calculations of the electronic structure and Tc for the 
case of single and multiple digital layers of GaAsMn and 
for the additional confinement effect by a quantum well 
of a Mn layer. We show that the behavior of the double 
layer as a function of the layer separation is essentially 
that of more than two layers. 

Our model is an extension to the case of planar het- 
erostructures of the Zener model Jox-bulk alloys of (III- 
V) ferromagnetic semiconductors. oll3 It is convenient to 
express the problem of the localized spins and the car- 
rier spins in an inhomogeneous systems in terms of the 
density functional theory.EHl It expresses the Curie Tem- 
perature in terms of the carrier spin susceptibility of the 
doped semiconductor, taking into account the spin-orbit 
interaction of the holes. Such an extension involves two 
technical refinements, compared to the bulk case. First, 
as confinement breaks the translational invariance along 
the growth axis z, the susceptibility becomes a non-local 
function x(z,z'). In order to determine Tc in a pla- 
nar heterostructure we have to solve an integral equation 
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whose kernel contains the non-local spin susceptibility. 
The derivation of the integral equation from the den- 
sity functional theory is given in section II. Second, as 
explained in section III, the calculation of the electronic 
structure of the holes in a planar heterostructure involves, 
the solution of the self-consistent Schrodinger equations 
The calculation of the subbands, for the in-planc_motion, 
is made using the k • p Luttinger HamiltonianJHj which 
includes the spin-orbit interaction. Our calculations ad- 
dress the regime in which both the Mn density, cm, and 
the carrier density, p, are high. This implies that several 
subbands are occupied. Therefore, our calculations in- 
clude both spin-orbit interaction-and multiple sub-bands, 
in contrast with previous work.e3~E3 

In section IV we apply the formalism to the case of a 
single digital layer of GaAsMn embedded in GaAs. From 
the calculated electronic structure for different values of 
the degree of compensation and the interdiffusion of the 
Mn, we find that the calculated Tq is an increasing func- 
tion the density of holes and a decreasing function of the 
interdiffusion of the magnetic impurities. In section V, 
we investigate the change in Tq as a function of the inter- 
layer distance between two layers in order to understand 
the observed behavior for multilayers of Mn for different 
degrees of compensation and interdiffusion. Our numeri- 
cal results for two layers, which are very similar to those 
for up to five layers, reproduce the main features of the 
experimental results. A remarkable result is found in sec- 
tion VI where the calculated Tc for a digital layer inside 
a quantum well is found to increase dramatically by the 
additional quantum well confinement, compared to just 
the confinement due to the Mn layer in section IV. In 
section VII we put our results in perspective and draw 
some conclusions. 



II. DENSITY FUNCTIONAL FORMULATION 
FOR MAGNETISM IN HETEROSTRUCTURES 

The observed effects of ferromagnetism in Mn-doped 
III-V semiconductors appear to be consistent with the 
microscopic mechanism of indirect Mn-Mn spin interac- 
tion mediated by the mobile holes via the exchange in- 
teraction between the hole and the magnetic moment of 
the localized d electrons of the Mn impurity. A com- 
mon model for the bulk system, which we shall adopt, 
consists of a quantum degenerate gas of fermions which 
interact, via a contact Heisenberg exchange interaction, 
with the local magnetic moments of Mn. The random 
array of Mn impurities is replaced by a homogeneous dis- 
tribution and the exchange interaction is reduced corre- 
spondingly. This procedure is dubbed the virtual crystal 
approximation.13 The exchange coupling between Mn and 
the holes is treated in the mean-field approximation.Ej 

The study of heterostructures of semiconductors 5- 
doped with manganese, calls for an extension of the vir- 
tual crystal and mean-field approximations to an inho- 
mogeneous distribution of Mn. We formulate the theory 



of magnetism of the inhomogeneous system in terms of 
the density functional theory, extendedptp include the 
spin densities and to finite temperatureEa The free en- 
ergy of the system of Mn spins and holes as a functional 
of the density and spin density distributions of the Mn 
spins and the hole carriers is separated into three contri- 
butions: 



F = F, 
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E h + E, 
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respectively, of Mn spins, the holes, and the Mn-hole in- 
teraction. The free energy of the Mn spin system alone 
is: 

F m [cm(t),M(t)] 
= [ d 3 rc M (r)[f (M(r)) + U M M 2 (r)} 



d 3 r / d 3 r'c M (r)c M (r')J MM (r,v')M{r)M(r% (2) 



The number density of the Mn is denoted by cjn(r) and 
M(r) is the spin expectation per Mn atom (in units of 
h). The noninteracting part of the free energy per Mn 
spin is given by 



/o(M) = k B T 



Mb-J^± 
\ sinh § 



(3) 



where kg is the Boltzmann constant, S = | is the spin 
of Mn and b — b(M) is related to the inverse of the usual 
Brillouin function, 



M = (S+ |) coth[(S* + \)b] - \ coth {\b) 



(4) 



The second term in the integral takes into account only 
the short-range Heisenberg exchange between Mn spins. 
Jm has the dimension-af energy. In GaAs it is believed to 
be antiferromagneticEj The third term accounts for the 
long-range dipole interaction Jmm{y, r') and is found to 
be negligible. 

The temperature range under study is sufficiently low 
compared with the Fermi temperature of the holes that 
the hole free energy will be taken as-jthe usual ground 
state energy functional Eh\p{v), S(r)]j23 where p(r) is the 
hole density and S(r) the spin density. The interaction 
between holes is included in the Hartree approximation in 
the calculation of the sub-bands. As the density of holes 
is very high, the exchange and correlation potential in 
the local-density approximation is of minor importance. 

The Mn-hole interaction term is given by 

%[p(r),S(r);c M (r),M(r)] 

d 3 r [ d 3 r'p(r)u(r - r')[c M (r') - c c (r')] 



(5) 



J J d 3 r c M (r)M(r) • S(r). 



The first term on the right side of the equation is the 
attractive potential provided by the Mn donors to the 
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holes with ju(r — r') being the Coulomb interaction. 
Experimentfl shows compensation of the acceptors, which 
are believed to be antisite impurities. The compensating 
impurity concentration c c (r) is taken into account. The 
second term is the hole spin interaction with the Mn spin, 
for which we use the simplest mean-field approximation. 
Functional terms beyond the mean field might be con- 
structed from theory such as Ref. [29|-|3l|. The hole-Mn 
spin interaction J has the dimension of energy- volume. 

In the mean-field and virtual-crystal approximations, 
the model describes the hole carriers interacting with an 
effective magnetic field produced by the localized Mn im- 
purities and vice versa. The variational result of the free 
energy functional with respect to both the magnetization 
of the Mn impurities and the magnetization of the hole 
carriers shows the interdependence of the two magneti- 
zations. Each is governed by the effective magnetic field 
generated by the other. They have to be determined 
selfconsistently. We report here only the work on the 
transition temperature. Theoretical finite magnetization 
studies are being carried out. Close to the Curie tem- 
perature, the magnetizations are small. The free energy 
is a quadratic functional of the two magnetizations. In 
a planar heterostructure with the growth axis along z, 
there is translational invariance along the x — y plane in 
the effective mass approximation so that quantities de- 
pend only on z. The free energy functional per unit area 
is then: 

F \p{z), S(z);c M (z), M(z)} = F\p(z),c M (z)] + 



dz c M {z) 



3k B T 



S(S + 1) 
J I dzc M {z)M a (z)S a (z) 



Jm 



M a {zf 



+ h dz / dz' S a (z)K a (z 7 z')S a {z'). 



(0) 



For simplicity, we investigate easy magnetization of 
M a (z) and S a (z) only along the growth axis a = z or in 
plane a = x. The first term F\p(z), cm{z)\ is the density 
functional for zero magnetization including the impurity 
potential for the holes. Charge neutrality determines the 
total number of holes: 



p(z)dz 



[cm(z) - c c (z)} dz. 



(7) 



The first quadratic term in M contains the inverse sus- 
ceptibility of the Mn spins. The last term is the magnetic 
energy of the holes, where Karif the inverse of the non- 
local hole spin susceptibility ;\I23: 



dy K a {z,y)x a {y,z') = S(z - z'). 



(8) 



Minimization of the energy functional (J6|) with respect 
to the magnetizations leads to two coupled equations for 
M a (z) and S a (z): 



c M (z) 



3k B T 



Jm 



+ JS a {z) \ =0 



(9) 



S(S + 1) 

dz' K a (z,z')S a {z') + Jc M {z)M a {z) = 0. (10) 



M a (z) = S a (z) = are always a solution, corresponding 
to the stable state only in the paramagnetic phase. The 
Curie temperature is the highest temperature at which 
Eqs. (|9|,[l0|) have non-zero solutions. Elimination of S a (z) 
from Eqs. (|Tj|,[l0|) by Eq. (||) leads to an integral equation 
for the Curie temperature Tq\ 



M a (z)c M (z) = 
c M {z)S(S+l)J 2 
3(k B T c + k B T M ) 



dzx a (z,z')c M {z')M a (z') (11) 



where k B T M = S(S + l)J M /3. Eq. (y), which relates 
Tq with the non-local spin susceptibility of the holes, 
Xa{z, z'), for a given planar heterostructure is the main 
result of this section. Eq. ( p"l] ) extends Eq. (10-12) of 
Rcf. to include multiple subband occupation and spin- 
orbit interaction. 



III. NON-LOCAL SUSCEPTIBILITY 

The calculation of the Curie Temperature for the pla- 
nar heterostructure involves the solution of the integral 
equation (jll]), whose kernel contains the non-local spin 
susceptibility. In this section we briefly describe our cal- 
culation of the electronic structure and the non-local spin 
susceptibility for a rather general planar heterostructure 
characterized by a profile of the Mn impurities cm{z), 
the profile of the compensating impurities c c (z) (anti- 
sites), and a band gap profile which creates a potential 
for the holes Vi(z). 

For a given density profile we solve the Poisson equa- 
tion and obtain the electrostatic potential V e i , which, to- 
gether with a band gap-potential Vi(z), defines the effec- 
tive mass Hamiltoniar&Ea for the envelope function of 
the holes: 



H, 



eff : 



H L (k 



ld_ 

' i dz 



V(z) 



(12) 



where k is the in-plane wave vector, V(z) — Vi(z) + 
V e i(z), Hl is the standard 4x4 Luttinger Hamiltonian, 
with parameters 71, 7 = 0.5(72 + 73), and [i = 0.5(73 — 
72). We adppt the cylindrical approximation about the 
growth axisB3, i.e., taking /x = in the Luttinger Hamil- 
tonian. This approximation has the advantage that 
does only depend on k = |k|. The physical properties 
involve angular integration over k about z so that the 
deviations from the cylindrical approximations are very 
smalltJ. 

The subbands ek,iy and the corresponding eigenstates, 
■0k, f (z), are given by the solution of the set of four cou- 
pled second order differential equations: 
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where 



ff eff Ck 



1^ 

i dz< 



V(z)6 r , 



IV. SINGLE DIGITAL LAYER 



F^(z) = £k ,^r (z) (13) 



,^3/2 



^1/2 ^-1, 1 



-,k,i/ 



(14) 



arc the four components of the wave function il^q V {z). 
Eq. (|l^) is solved with the mini-band k • p methodEl We 
first solve the |k| = case, in which the equations are 
decoupled into two ordinary Schrodinger equations cor- 
responding to the light and the heavy holes. These have 
rift and ni bound states, evaluated by transforming the 
one dimensional Schrodinger equation into a tridiagoual 
matrix eigenvalue problem which is solved numerically^! 
The |k| = solutions form a basis set with the 2nh 
states i/) u hO-, 0, 0, 0), and ip u h{0, 0, 0, 1) and the 2m states 
^1/2(0,1,0,0) and Vfi(0> 0, 1, 0). The finite |k| eigenergies 
£k,i/ and eigenstates ipk lV (z) are obtained in terms of the 
basis set as the solutions of Eq. ( |l3| ) as the N x N secular 
determinant problem, where N = 2(nh + ni). The hole 
density is given by: 



p( z ) = J2 J 7^4/Ki/)^k,i/(^ 



(15) 



where /(e) is the Fermi-Dirac occupation function, the 
Fermi level being fixed so that the charge neutrality con- 
dition is met. The hole density p(z) and the potential 
V(z) are determined by iteration to self-consistency. 
The nonlocal spin susceptibility is then given by: 



X, 



d 2 k 



d/( e k,i/) a . , Q - 
9ek ) y ' °\l,v,v\ z ) 



i/'^i/ 



£!/ (k)-ei//(k) 



, (16) 



where a = (x,y,z) and the spin matrix elements are 
given by: 



SS v ,Jz) = (^y(z)\S a \i; k<v (z)), 



(17) 



where the angular brackets denote the expectation value 
over the spin degrees of freedom of the hole states ipk, v (z). 
In the absence of spin-orbit interaction, the spin ma- 
trix elements would be independent of the in-plane mo- 
mentum k. Moreover, due to the interplay between the 
spin-orbit interaction and confinement, the nonlocal spin 
susceptibility takes different values for in-plane and off- 
plane orientation. If the spin-orbit interaction was the 
only source of anisotropy, our calculation could deter- 
mine the easy axis. Other sources of anisotropy, like 
shape anisotropy, are not considered in our calculation. 
In the calculations reported in this paper, we assume 
an ini-»lane magnetization, guided by the experimental 
resulted. 



In this section we present the results of our calcula- 
tions of the electronic structure and the Curie tempera- 
ture, Tc, for GaAs doped with a single digital Mn layer 
of (Gao.5Mno.5)As, as in the experiments. In an ideal 
(5-doping, the Mn atoms occupy a single atomic plane. 
In the real system, Mn atoms undergo interdiffusion to 
occupy several layers. We assume that the compensating 
impurities are closely associated with the Mn atoms and 
assume their distributions to have the same shape: 



cm(z) _ c c (z) 



1 



OLM 



^ A 



e -(™/^) 2 S( z -na), (18) 



where a c and aju lead, respectively, to the total concen- 
trations of Mn, cm and of the compensating impurities, 
c c so that the density of holes is p = cm — c c . Hence, for 
a given Mn concentration cm, a single layer is character- 
ized by (A,p) 

In the limit of A = 0, we recover the ideal 5-doping 
case, cm(z) — cqS(z). Then, Eq. ([ll]) can be solved ana- 
lytically: 



k B T c = \c M S{S + 1) J 2 X (0, 0) - k B T M . 



(19) 



The dynamics of the holes is contained in %(0, 0). Wc 
see that Tc does not depend on the sign of the hole Mn 
exchange interaction, J. On the other hand, — ^b^m, 
proportional to the direct Mn-Mn interaction, if antifer- 
romagnetic (Jm > 0), decreases the Curie temperature, 
as expected. In the case of a random alloy of GaAsMn, 
Jm is found negligibleE-3 because the distance between the 
Mn impurities is rather high. In contrast, the in-plane 
average distance between the Mn is much shorter in the 
digital heterostructure. However, an accurate -value for 
both J and Jm is not known. Hereafter, we setfl Jm = 
and J = 150 mcV-nm 3 . 
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FIG. 1. Left panel: selfconsistent potential for a single dig- 
ital layer with A = 0.5 nm and p = 1.3 x 10 13 cm -2 , includ- 



ing the HH and LH levels. Right panel: 
dashed line is the Fermi level. 



hole subbands. The 



In the left panel of Fig. [I] we plot the selfconsistent po- 
tential for the holes corresponding to a single digital layer 



4 



with A = 0.5 nm and p = 1.3 x 10 13 cm" 2 , together with 
the energy levels for the light and the heavy holes. In the 
right panel we plot the subbands for the in-plane motion 
of the holes. The dashed line indicates the Fermi level. 
For Gao.5Mno.5As, we have cm = 3.13 x 10 14 cm -2 . Even 
for a density of holes at only 4.1% of the Mn concentra- 
tion, 3 subbands are occupied with holes. For this set 
of parameters the obtained T c is 35-Kpthe experimental 
value obtained by R. Kawakami et al.tB 

The spin-orbit interaction causes both the anticross- 
ing and the non patabolic shape of the hole subbands. 
As in the bulk casc,ll3 the spin-orbit effect also reduces 
significantly the effective magnetic coupling between Mn 
spins. Therefore, it is important to include spin-orbit 
in the theory of ferromagnetism in planar heterostruc- 
tures, an ingredient missing in previous papers for quan- 
tum wells. tB Ej 
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FIG. 2. Tc as a function of density of holes, p, for A =0, 
0.5, 1.0, 1.5 and 2.0 nm (from top to bottom). Inset: j3 as a 
function of A for the fit T c (p, A) oc p /3(A) . 
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FIG. 3. Tc as a function of the interdifusion parameter, 
A, for p — 1.5 (highest), 0.6, 0.3, 0.15 (lowest) in units of 
10 14 cm" 2 . Inset: Lines of Tc= 35 K (crosses) and Tc =60 K 
(stars), in the (p, A) plane (with J = 150 meV-nm 3 ). 

In Fig. H we show Tc for a single layer, as a function 
of the density of holes, for different values of the inter- 
diffusion parameter, A, and for a fixed value of the Mn 
concentration cm — 3.13 x 10 14 cm" 2 . We take the mag- 
netization to be in the plane of the layers in line with 



the experimental result.ta However, preliminary calcula- 
tions show that the model predicts an off-plane easy axis 
for the single digital layer, in the idealized A = case. 
For each point we have calculated the electronic struc- 
ture selfconsistently and solved Eq. JTl| ) to obtain Tc- 
The general trend is that Tc is an increasing function 
of the density of holes. The lines are the best fit using 
T c oc pP. In the inset we plot j3 as a function of A. In 
a two-dimensional system with parabolic bands, no spin- 
orbit interaction, and no Coulomb interaction, we would 
have obtained j3 = 0. Remarkably (3 in the siaale layer 
is even larger than the value for bulk, /3 = 1/3.E-3 

In Fig. |^ we plot Tc for the same set of single layers, 
as a function of A, for different densities of holes. The 
model shows that interdiffusion reduces Tc- For a fixed 
value of J, there is a line in the plane (p, A) which gives 
the same Tc- In the inset of Fig. || we plot that line 
for both T c = 35 K and T c = 60 K for J = 150 mcV- 
nm 3 . The first (35K) corresponds Xo the Curie temper- 
ature reported by Kawakami et alUB to yield an idea of 
what model parameters could describe the experimental 
conditions. The second (60K) has been obtained for the 
same kind of— heterostructures grown at slightly higher 
temperature] 2 .!! 

V. DOUBLE DIGITAL LAYER 

In this section we report on our calculations of the elec- 
tronic structure and Tc for two identical digital layers, 
separated by N monolayers of GaAs, so that the inter- 
layer distance is d = N x 0.2825 nm. Both layers are 
described by Eq. (p8[). We choose a point in the (A,p) 
parameter space, so that, for very large d, the calculated 
Tc is close to the experimental value of 35 K. Then we 
calculate Tc(d) for smaller values of d. 
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FIG. 4. Round dots: Curie Temperature as a function of 
the interlayer distance for a double layer system. In the left 
panel, A = 0.5 nm and p — 1.3 x 10 13 cm -2 . In the right 
A = 1.5 nm and p = 3 x 10 13 cm" 2 . Square |dpts: experi- 
mental Curie temperature for a multilayer caseO 

The results are shown in Fig. ^ for two cases: (A = 
0.5 nm, p = 1.3 X 10 13 cm" 2 ) (left panel) and (A = 
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1.5 nm, p = 3 x 10 13 cm -2 ) (right panel). The theo- 
retical results obtained with-the first case give a better 
fit to the experimental datall3 than those obtained with 
the second. In Fig. || we plot the corresponding density 
profiles for 3 different interlayer distances to represent 3 
regions of separation dependence in the Curie tempera- 
ture. At short layer separations (upper panels of Fig. [^), 
both layers of the hole and of the Mn distribution overlap 
and T c depends strongly on the separation d. At inter- 
mediate separations (left middle panel of Fig. ||) , the two 
layers of Mn do not overlap but the hole distribution still 
does. As a result, the layers are coupled and Tc is weakly 
dependent on the interlayer distance. A further increase 
of d leads to the uncoupled regime where Tc reaches the 
single layer value (lower panels of Fig. ||) . 
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FIG. 5. Mn (shaded regions) and hole density (lines) pro- 
files for a double layer. Left panel: p — 1.3 x 10 13 cm -2 per 
layer and A = 0.5 nm. Right panel: p — 3 x 10 13 cm -2 per 
layer and A = 1.5 nm. From top to bottom, the interlayer 
distance is 10, 20 and 40 monolayers. 

Both the calculated and the measured Tq decrease as 
d increases and they reach a stationary value at large d. 
Similar behavior is obtained for several values of (p, A). 
We have also calculated Tc for 3, 4, and 5 delta layers 
and the separation dependence is similar. The steep de- 
cline of Tc with separation stops at interlayer distances 
higher than about 10 monolayers, the relevant experi- 
mental region. In all these cases, the theory seems to 
slightly underestimate the coupling between the layers at 
intermediate distances. This might indicate that the den- 
sity of itinerant carriers in between the magnetic layers 
might be higher in the experiment than in our calcula- 
tions. Further theoretical and experimental work might 
clarify this point. 

VI. DIGITAL LAYER INSIDE A QUANTUM 
WELL 

In this section we study the electronic structure and 
the Tc of a single digital ferromagnetic layer inside a 
quantum well, using the formalism of sections II and III. 
The model predicts that the confinement effect increases 



Tc up to a factor of 3 compared with the unconfined 
single layer. The system consists of Ga 2; AsAl 1 _ :I , barri- 
ers containing a GaAs well with a single digital layer of 
Gao.5AsMno.5 in the middle. This structure would be the 
ferromagnetic -analog of the Be (S-doped GaAs/GaAsAl 
quantum well.L3 
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FIG. 6. Hole density profiles for different values of the bar- 
rier potential, V, for a Mn layer (vertical lines) inside a quan- 
tum well. The total density of holes is p = 1.3 x 10 13 cm -2 . 
For the sake of clarity, the Mn distribution has been divided 
by 20. 

We model the interface GaAs/Ga^AsAli-a, as a bar- 
rier potential of height V = 550 • x meV. The distri- 
bution of Mn is that of section IV, with A = 0.5nm, 
p = 1.3 x 10 13 cm~ 2 . In Fig. || we show the distribu- 
tion of holes of a 10 monolayers wide quantum well for 
different values of the aluminum content. As the barrier 
potential increases the hole distribution overlaps more 
the Mn layer, increasing the Tc ■ In Fig. [?| we plot 
Tc as a function of the barrier height (proportional to 
the Al content of the barrier), for a quantum well width 
of both 10 and 20 monolayers (2.8 and 5.7 nm) with a 
Mn layer in the middle of the well, with (A = 0.5 nm, 
p = 1.3 x 10 13 cm -2 ). For the narrower well the enhance- 
ment factor can be as large as 2.9 for V = 300 meV which 
corresponds to 54% of Aluminum in the barriers. For a 
wider quantum well the effect is smaller. 

The increase of Tc is due to two factors. The first 
is the increase in the overlap between the Mn and 
the hole distributions (see Fig. |). In Fig. 0b we 
plot the overlap of the Mn and the hole distributions, 
/ dzy / p(z)cM(z) I \jpcu , as a function of the barrier 
height, V. The increase of the overlap due to confine- 
ment is larger for the narrower well. The second factor 
is the increase of the density of states at the Fermi Level 
(DOS). In figure ^c we plot the DOS as a function of V, 
normalized by the DOS at the Fermi level for the case 
V = 0. At V = the Fermi Level is close to the bottom 
of the second heavy hole band so that there 3 bands are 
occupied (see figure 1) . The effect of the barrier poten- 
tial is to increase the energy level spacing so that as V 
increases, the Fermi level goes below the second Heavy 
hole band and moves towards the bottom of the first light 
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hole band. For small values of V the DOS at the Fermi 
level decreases slightly, increasing up to 40% for higher 
values of V. 
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FIG. 7. (a) Curie Temperature of a single digital layer in a 
quantum well as a function of the barrier height, V, for two 
different values of the well width, Lqw- (b) Overlap of the 
hole and the Mn distributions as a function of V. (c) Density 
of States (DOS) at the Fermi level (in units of the DOS for 
V = 0) 

For V > 350 meV and Lqw =10 nm, the Fermi Level 
gets close to the bottom of the light hole band where 
there is a dramatic increase of the DOS at the Fermi 
Level. This leads to an even larger increase of the pre- 
dicted Critical Temperature. Work is in progress to check 
if this result (not shown in the figures) remains when the 
finite temperature of the fermions is taken into account. 

We have also checked that, for a 10 monolayer quantum 
well with a density of holes of p — 2.5 x 10 13 cm -2 , the 
relative increase of Tq for 1^=300 meV, is a factor of 1.8, 
i.e., smaller than the enhancement for p = 1.3 x 10 13 

_9 

cm ~ 



VII. DISCUSSION AND CONCLUSIONS 

The work presented in this paper is based on a model 
of effective-mass, virtual- crystal and mean-field approxi- 
mations. An important feature of the theory is that Tc 
scalesjpijth the square of the exchange coupling constant, 
jJiSBB We have used a value of J = 150 meV -nm from 
Ref. HJI^jJrlowever, a value 3 times smaller has also been 
reported.LS This would change Tc by a factor of nine. In 
turn, this could be compensated in the theory by increas- 
ing the density of holes, including the split-off band in the 
calculation,E3 a larger effectrpj mass associated with the 
motion in an impurity band,E3 or a direct ferromagnetic 
coupling Jm between the Mn atoms. 

These caveats indicate the current difficulties with the 
theoretical prediction of the absolute values of the Curie 
temperature. Our model is less ambitious and is used to 
provide changes in the behavior of the Curie temperature 
as a function of parameters which characterize the pla- 
nar heterostructure. For instance, we have presented an 



analysis of Tc for a single digital layer as a function of 
the density of compensating impurities and the spread of 
the Mn atoms due to interdiffusion. Our results indicate 
that Tq increases with the density of holes and with less 
interdiffusion (smaller A). The fact that Tc depends on 
density in a quasi-two dimensional system, contrary to a 
naive calculation with parabolic bands, gives some theg 
oretical support to the experiment by H. Ohno et al.tH 
in which they observed a change of Tc as the density of 
holes changes in a InAsMn quantum well in a field-effect 
transistor. Our theory could be also applied to model 
ferromagnetism observed in a quantum well of p-doped 
CdTeMnB 

In section V we have presented our calculations for the 
double layer system together with the experimental re- 
sults for multilayers. The qualitative agreement is good, 
but the theory seems to underestimate the Tc at inter- 
mediate interlayer distances, i.e., the coupling between 
the magnetic layers is larger in the experiment than in 
the theory. Further work on this problem, both in exper- 
imental characterization and in the improvement of the 
theory, might shed some light on the microscopic origin 
of ferromagnetism in this kind of systems. 

In section VI we have studied the confinement effects 
on the hole carriers which mediate the Mn-Mn magnetic 
coupling. This leads to a prediction of an increase of 
Tc by as much as a factor of almost 3, when a single 
digital layer is grown in a quantum well structure. The 
capability to investigate the systemic changes on ferro- 
magnetism and to predict observable effects is a strong 
point of the effective-mass mean-field theory. 

In conclusion, we have presented a theoretical frame- 
work to calculate the electronic structure and the critical 
temperature of heterostructures of III-V ferromagnetic 
semiconductors. This work is an extension of the 3D 
case,t3 in which the relevance of the spin-orbit interac- 
tion has been pointed out. The main features of the for- 
malism, absent in previous papers on heterostructures, 
are the inclusion of several subbands, necessary because 
of the high density of holes in the system, and the in- 
clusion of the spin-orbit interaction, important because 
it changes both the magnetic coupling and the shape of 
the bands. We have presented calculations of the digital 
magnetic heterostructures, providing a qualitative under- 
standing of the experimental values of Tc, and we have 
predicted that Tc for a single digital layer can increase 
by a factor of 2 when embedded in a quantum well. 

We wish to thank Drs. D.D. Awschalom, R. Kawakami, 
and A. Gossard for stimulating discussions and Dr E. 
Gwinn for suggesting the calculation of section VII. 
We acknowledge Spanish Ministry of Education for a 
post-doctoral fellowship and support by DARPA/ONR 
N0014-99-1-109 and NSF DMR 0099572. 
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